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BSTRACT  (Continue  on  teveree  side  II  necessary  end  Identify  hy  block  number) 

An  ultrasonic  viscometer,  capable  of  measuring  small  liquid  volumes  such  as 
drops,  was  developed  by  utilizing  the  short  effective  path  length  of  shear 
waves  in  liquids.  The  viscometer  consists  of  an  unplated  quartz  crystal 
excited  to  shear  vibrations  with  the  liquid  drop  to  be  tested  placed  on  its 
surface.  The  vibrations  of  the  natural  modes  of  the  quartz  plate  was  repre- 
sented by  a mass-point  compliant  equation,  and  the  viscosity  equation  was 
derived  from  the  knowledge  of  the  energy  loss. 

Suitable  electronics  was  built  to  drive  the  crystal  and  to  record  the  decayin 
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20.  ABSTRACT  (continued) 

-A  vibrations.  A crystal  orientation  analysis  insured  a dominant 

thickness-shear  mode  not  coupled  to  flexure  and  face-shear  motion. 
The  amount  of  sound  radiated  was  determined  from  the  change  In  the 
quality  factor  when  the  crystal  was  vibrating  alone  In  air  and 
».der  vacum.  ^ 

Heasuren^nts  were  made  on  a variety  of  normal  liquids  and  bonding 
agents  over  the  frequency  range  from  200  to  1000  kHz,  and  the 
dynamic  viscosity  results  scattered  about  the  static  value.  A 
number  of  body  fluids  of  man  and  domesticated  animals  were  examined. 

The  response  of  the  shear  viscometer  to  changes  in  the  rheological 
properties  of  bovine  red  b.'ood  cells  was  also  measured.  That  Is, 
the  ability  of  red  cells  to  deform  was  changed  by  hematological 
methods  and  the  voscoslty  was  measured  before  and  after.  It  was 
found  that  the  viscosity  Increased  with  the  addition  of  sodium 
oleate,  acetaldehyde,  and  hematocrit. 
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CHAPTER  I 


INTRODUCTION 

At  the  present  time,  e3q>erimental  information  dealing  with  the 
measurement  of  the  viscosity  of  small  liquid  volumes,  such  as  drops,  is 
not  available  between  150  kHz  and  5.0  MHz  (1).  Probably  the  first  work 
done  in  this  direction  was  performed  by  Skudrzyk  and  Dussik  in  the 
hospital  at  Ischl  in  1947.  A thin  layer  of  the  fluid  found  in  the 
ventricle  cavities  of  the  human  brain  was  placed  over  the  surface  of  a 
quartz  crystal  that  was  excited  to  dilatational  vibrations  at  its  first 
three  harmonic  frequencies  of  100  kHz.  The  fluid  increased  the  damping 
of  the  crystal  considerably  and  it  seemed  as  if  the  large  albumin 
chains  in  the  fluid  were  responsible  for  the  relaxation  phenomenon  in 
the  frequency  range  from  100  to  500  kHz. 

In  1968,  Ghering  (2)  continued  this  work  ar»d  investigated  the 
energy  losses  in  free  and  bonding  liquids.  The  experimental  data  for 
films  was  puzzling;  the  magnitude  of  the  energy  losses  were  much 
greater  than  could  be  explained  on  the  basis  of  longitudinal  vibrations. 
The  discrepancy  was  finally  explained  by  the  presence  of  shear  motion 
in  the  quartz  disc  vibrator.  Such  motion  usually  leads  to  greater 
losses  in  free  liquids.  However,  the  energy  dissipation  for  loaded 
films  sandwiched  between  two  transducers  was  still  too  great  even  when 
shear  motion  was  accounted  for. 

Hosier  (3)  continued  and  extended  the  methodology  using  a shear 
vibrator.  Hc«fever,  he  encountered  numerous  difticuities  that 
prevented  the  author  from  obtaining  more  detailed  results.  The  range 


of  difficulties  can  best  be  illustrated  by  the  following  quote  (4) : 

. the  values  of  viscosity  measured  for  water  scattered  between 
0.5  centipcise  to  450  centipoises.  Similar  results  were  obtained  for 
glycerin  where  the  viscosity  ranged  from  400  to  1100  centipoises  where- 
ais  the  static  viscosity  is  approximately  600  centipoise." 

1. 1 Statement  of  the  Problem 

It  is  through  these  studies  that  a method  to  measure  the  viscosity 
of  free  liquids  became  possible.  The  task  of  this  study  was  to 
eliminate  the  experimental  errors  of  the  previous  authors  and  develop 
a high  precision,  inexpensive  viscometer  capable  of  measuring  liquids 
of  small  volume,  such  as  drops.  This  study  utilizes  the  short  effec- 
tive path  length  of  shear  waves  in  liquids  to  examine  drops.  In 
liquids,  shear  waves  are  strongly  attenuated  and  generate  a boundary 
layer  of  motion  in  much  the  same  way  as  electromagnetic  waves  generate 
the  skin  effect  in  metals,  where  an  effective  resistance  similar  to 
boundary  viscosity  exists-  By  analogy. 


Y = 


tujpj 


a 


(1.1) 


where  Y is  the  skin  depth  defined  as  the  propagation  distance  required 
for  the  amplitude  of  a shear  wave  to  fall  to  1/e  of  its  initial  value. 
In  Equation  (1.1) 

y = viscosity  in  poises, 

u)  = angular  frequency  in  rad/sec, 

p = density  in  gm/cm^,  and 

a = absorption  coefficient  in  cm 
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Table  I illustrates  the  variation  of  the  propagation  distance  y,  the 

shear  wavelength  X auid  the  shear  velocity  v , , for  water  over  the 
sh  sh 

frequency  range  from  0.2  to  1.0, MHz.  Typically,  in  this  frequency  range, 

, . '•  a i' 

a shear  wave  is  extinct  after  traveling  a distance  of  1.0  micron. 

A method  for  measuring  the  viscosity  of  drops  of  liquids  at  high 
frequencies  would  have  numerous  practical  applications.  The  stemdard 
method  of  measuring  viscosity  at  the  low  physiological  strain  rates  is 
to  use  a Brookfield  Couette  Viscometer  which  requires  500  milliliters 
of  liquid.  Equally  significant,  a thin  layer  of  the  test  liquid  does 
not  slip  along  the  rotating  cylindrical  wall  and  could  possibly  conceal 
the  true  viscosity  if  a mixture  were  measured. 

At  the  present  time,  no  viscometer  is  capable  of  measuring  the 
viscosity  of  body  fluids,  such  as  those  found  in  the  joints,  spinal 
column,  and  the  ventricular  cavities  of  the  brain.  A simple,  direct 
measurement  of  their  viscosity  could  be  indicative  of  their  state  of 
health  and  diminish  the  need  for  existing  operative  techniques. 

A droplet  viscometer  would  be  a new  tool  for  examining  other 
liquids  of  special  interest  to  researchers  in  the  medical  fields.  No 
viscometer  is  presently  capable  cf  measuring  uncoaguiated  whole  blood 
in  less  than  four  minutes.  The  shear  viscometer  makes  it  possible  tc 
measure  the  viscosity  of  fresh  blood  from  a finger  prick  in  less  than 
one  minute. 

Medical  schools  are  also  interested  in  the  effect  that  a change  in 
the  ability  of  red  blood  cells  to  deform  has  upon  the  viscosity  of 
blood.  The  fluid  nature  of  red  blood  cells  and  their  incomplete 
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Table  I.  Shear  parameters  of  water 


F(kHz) 

Y(cm) 

a (cm' 

■^) 

(cm) 

fern] 

'^sh  1 s J 

235 

1.16 

X 

10”** 

8.59 

X 

10^ 

7.30 

X 

10-“ 

172 

283 

1.06 

X 

10-“ 

9.43 

X 

10^ 

6. 66 

X 

10-“ 

188 

402 

0.88 

X 

io"“ 

11.23 

X 

10^ 

5.58 

X 

10-“ 

224 

503 

0.80 

X 

10"“ 

12.57 

X 

10^ 

5.00 

X 

10-“ 

251 

1000 

0.56 

X 

10"  “ 

17.70 

X 

10^ 

3.54 

X 

10"  “ 

354 

filling  allows  normal  red  cells  to  be  deformed  in  a variety  of  shapes, 
ranging  from  a disk  to  a spheroid.  Consequently,  a deformed  red  cell 
will  respond  differently  to  a shear  wave  through  the  flow  of  the 
internal  fluid  of  the  cell  aind  the  bending  of  the  outer  membreme.  This 
leads  to  changes  in  viscosity. 

The  study  of  erythrocytes  is  very  pertinent  as  an  application  of 
the  shear  viscometer,  since  several  authors  have  found  the  response  of 
red  blood  cells  to  shear  forces  to  be  qualitatively  similar  to  those  of 
fluid  drops  in  suspension  (5). 

1. 2 Quartz  Shear  Viscometer 

The  shear  viscometer  consists  of  em  unplated  quartz  crystal 
operating  at  the  fundamental  frequency  and  excited  to  shear  vibrations 
by  ultrasonic  pulses.  The  unobstructed  top  surface  is  in  contact  with 
the  liquid  drop,  and  the  bottom  surface  is  returned  to  the  hot  side  of 
the  generator.  After  the  power  is  shut  off,  the  crystal  becomes  a 
receiver  and  converts  the  decaying  vibrations  into  an  electrical  signal. 
These  vibrations  are  amplified  by  a high  gain,  low  leakage  pulse 
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eunplifier  and  displayed  on  a storage  oscilloscope.  The  decay  time  of 
the  vibration  is  determined  from  the  oscilloscope  curve. 

Because  of  the  great  dairying  generated  by  the  fluid,  it  was  not 
necessary  to  design  precision  crystal  holders  to  reduce  interaction 
losses  into  the  suspension.  The  quartz  was  simply  placed  on  a flat 
metal  base.  The  air  cushion  was  sufficient  to  isolate  its  vibration 
from  the  base.  The  second  electrode  is  a stiff  metal  plate  located 
about  5.0  nun  ^d^ove  the  quartz  surface,  but  not  exactly  parallel  to  it. 
Because  of  this  separation  wd  inclination,  st£mding  waves  in  the  air 
gap  were  of  no  consequence  to  the  vibration  of  the  face  of  the  crystal. 
On  the  other  hand,  this  type  of  crystal  holder  confined  sound  radiation 
to  the  air  spaces  around  the  crystal  amd  led  to  a quality  factor  of  the 
order  of  10®  in  air. 

Because  of  the  finite  size  of  the  quartz  crystal,  bending,  torsion, 
shear  and  dilatational  displacements  are  closely  coupled  and  a quartz 
vibrator  usually  has  many  resonemces  close  to  the  main  mode  of  vibra- 
tion (see  Appendix  A) . Damping  may  then  affect  the  main  resonance 
differently  from  its  neighbors  and  the  vibration  that  appears  ;n  the 
decay  may  be  very  different  from  the  one  that  had  been  excited  by  the 
generator. 

The  main  problem  here  consists  of  finding  a quartz  crystal  that 
has  its  resonaince  frequencies  well  separated  from  spurious  frequencies 
so  that  the  mode  that  dominates  the  decay  is  also  the  mode  that 
generates  the  mam  resoneince  peeik  of  the  undamped  crystal.  It  turns 
out  that  AT  cut  quartz  bars  made  of  certain  dimensions  isee  Appendix  B) 
lead  to  very  satisfactory  frequency  response  curves. 


the  ciystal,  whereas  the  small  amplitudes  of  vibrations  are  the  limit- 
ing factor  in  the  upper  frequency.  The  viscometer  is  also  restricted 
as  to  the  aunount  of  fluid  that  can  be  placed  on  its  surface.  Normally, 
a sample  thickness  of  2.0  to  5.0  mm  is  used.  Increasing  the  liquid 
thickness  beyond  this  ramge  serves  no  purpose,  since  all  the  shear 
losses  occur  at  the  bottom  of  the  layer  within  a distance  of  1.0 
micron.  A thicker  liquid  will  not  increase  the  viscosity,  since  a 5.0 

mm  thickness  already  acts  like  w infinite  medium. 

t 


CHAPTER  II 

THEORY  OF  THE  SHEAR  VIBRATOR 

The  vibration  of  an  infinitely  extended  quartz  surface  (in  the  x 
direction)  vibrating  in  the  thickness-shear  mode  is  described  by  the 
standing  waves 

5 = C®,  cos  k y , (2.1) 

V J 

where  is  the  constant  ^^nplitude  and  k = mr/2b.  If  the  dimension 
along  y is  very  great,  standing  shear  waves  along  the  y direction  are 
also  solutions.  These  inodes  are  known  as  the  zero-order  uncoupled 
modes  for  infinitely  extended  surfaces.  For  vibrations  of  finite 
quartz  plates,  one  shear  mode  cannot  satisfy  the  boundary  conditions 
at  the  side  surfaces,  emd  other  modes  of  vibration  exist  and  are 
coupled  to  each  other,  each  one  generating  the  other  type  of  motion. 
The  solution  of  a long  quartz  bar  is  always  approximated  by  a sum  of 
such  zero-order  mode  functions  (6)  with  different  values  of  k, 

n 

C = I cos  k^  y . (2.2) 

The  superposition  of  the  modes  are  spread  out  in  a finite  wave  number 
range  k^,  . . . . . ' leads  to  a beat-like  amplitude  distribu- 

tion in  space, 

K = K'  (x)  cos  k y , (2.3) 

0 

where  ^'(x)  is  a slowly  varying  function  of  x,  with  a maximum  value  in 
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the  center  and  smaller  values  at  the  ends.  However,  the  exact  velocity 
distribution  'x)  for  quartz  surfaces  of  finite  extent  are  unknown  and 
must  be  found  experimentally. 


The  vibration  of  the  natural  modes  of  a quartz  plate  can  also  be 
represented  by  eui  equation  similar  to  that  of  a point  mass  vibrating 
with  a resistance  and  compliance  (7),- 


" ''v'v  * 7 

K 

V 


(2.4) 


it  It  it 

The  quantities  are  the  mode  parameters  auid  represent  the 

effective  mass,  resistance,  and  compliance  for  a particular  natural 
function  5^,  referred  to  the  observation  point  A and  the  total  driving 
force  F.  The  mode  parauneters  are  given  by 

* M (A) 

* ^ 

R (A)  “ -7  TlM^  (A)  (2.5) 

V (Jj  V 

and 


K^(A) 


% %(A) 


f 


where  H (A)  is  the  fraction  of  the  total  driving  force  that  is  availadDle 
V 

th 

for  exciting  the  v mode,  and  is  known  as  the  excitation  constant. 


If  the  driving  point  coincides  with  the  point  of  observation, 

= 1 regardless  of  the  complexity  of  the  system  (8) , then  the  effec- 
tive mode  parameters  become 
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<r^> 

V 

Mv M 


e. 


0) 


(A) 


R = 
V 


— riM 
0)  V 


emd 


(2.6) 


U)'-  M 
V V 

where  are  the  natural  functions  and  the  corresponding  natural 
frequencies. 


In  our  case,  the  force  is  removed  and  the  free  vibrations  decay 
with  time  because  of  internal  friction.  Experiments  show  that  the  decay 
is  exponential  when  damping  is  not  excessively  large  (9).  Differen- 
tiating and  rearranging  terms,  letting  the  natural  mode  now  be 
expressed  by  C, 


dt 


1 2 • d • n • 

L-  ^ + 26  ^ C + cdH  = 0 , 

2 dt  V 


where 


6 = — — = decay  constant 
2M  . 


and 


(2.7) 


resonance  frequency 


A solution  is  given  by 
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and 


where 


(2.8) 


* 

(ii 

0 


and  C has  a periodic  amplitude  decay  given  by 


(2.9) 


0)  n 

2_  f. 

• -6t  2 ^ 

^0  ® “ ® 


(2,10) 


where  6 is  the  decay  constant  of  the  envelope  of  the  vibrations  emd 
* 

oJj  is  the  effective  resonant  frequency  of  the  point-mass  element. 


Assuming  that  the  maximum  displacement  is  equal  to  the  displace- 
ment amplitude,  we  have 


K ==  K.  (-<5)e 


■6t 


(2.11) 


The  maximum  elastic  energy  of  the  vibration, 


E - i 


(2.12) 


is  then  given  by 


„ 1 w x2r2  -26t 

E - 5 5,  ^ 


dE 

The  rate  of  energy  lost  per  cycle  ~ is 

at 


(2.13) 


■7^ 
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= - 26e 

dt 

emd  the  decay  then  follows  from 
dE 


2E 


(2.14) 


(2.15) 


The  decay  constemt  is  therefore  a function  of  the  energy  loss  amd  the 
maucimum  vibrational  energy  of  the  crystal.  Because  of  the  strong 
coupling  of  the  dilatational  amd  shear  modes  due  to  the  boundary 
conditions  at  the  side  surfaces,  E represents  the  total  elastic  energy, 
consisting  of  shear  and  dilatational  components,  i.e.. 


E 


X V 2 

2 sh  sh 


V 2 
M V 
dll  dil 


(2.16) 


where  V , amd  V.,,  are  the  velocity  amplitudes  for  shear  amd  dilatation, 
sn  all 


Assuming  , 

sn  axi. 


E * i mV.  (1  + A*) 

2 V sh 

where 


(2.17) 


A = 


V 


dil 


sh 


(2.18) 


Let  us  consider  the  semi-infinite  fluid  layer  (y  Z 0)  which  is 
excited  to  shear  vibrations  by  the  plane  crystal  surface  at  y = 0.  The 
surface  of  the  crystal  moves  in  its  plane  sinusoidally  with  am  ampli- 
tude V . = V amd  a frequency  (*).  Using  Newton's  second  law  for  the 
sh  0 

resultant  stress  on  a layer  of  fluid  yields  the  Stokes  wave  equation. 


The  velocity  field  for  periodic  waves  is  described  by 


where 
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dv  ,,  3^v 
3y^ 


(2.19) 


-k^y  j(cat-kj|y) 

V = V e e 
0 


(2.20) 


UyJ 


(2.21) 


The  reaction  force  exerted  by  the  fluid  due  to  this  motion  is  given  by, 


^ y*0  “ 3yl„_ 


y«0 


= yk  V a (cos  wt  - sin  cot)  , (2.22) 

0 0 L 

and  the  power  dissipated  cthat  is,  the  work  performed  by  the  crystal 
surface  on  the  viscous  film)  is 


i? 

3t 


rT 


y=0 


= lim 


uk  V (cos  cot  - sin  cot)  V cos  cot  dt 
0 0 0 


^ yk  0 

2 0 L 0 


(2.23) 


The  decay  constant  thus  becomes 


dt 

2E 


(2.24) 


1 yk  0 

2 0 L 0 


2 i mV(1  + A2) 

2 V 0 


i 


ppp 


and 
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6 = 


yk  a 
2 0 L 

M®(1  + A^) 


(2.25) 


Since  the  increase  of  the  decay  constant  for  a crystal  loaded  by  a drop 
of  fluid  is  linearly  proportional  to  the  increased  shear  loss, 


I 

1 lUlD  1^ 

± u rr.  a 

2 ^ UyJ  L 
M^(l  + A^) 


(2.26) 


and 


2w  (An 


(1  + A*)^ 


(2.27) 


'■^e  viscosity  of  a small  volume  of  liquid  depends  on  the  chemge  in 

2 


system  loss  factor  An  = 


oift**  - t I 

Vo  c' 


- , the  liquid  contact  area  0^,  the 


crystal  frequency  u),  the  liquid  density  p , the  correction  factor 

L 

(1  + A^)^  that  accounts  for  the  presence  of  dilatational  motion,  and 
the  crystal  mode  mass  M®.  All  of  these  parameters,  except  the  liquid 


density,  must  be  determined  experimentally. 


ilSP31!PP*piiiP! 
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CHAPTER  III 

E.'TPERIMENTAL  EQUIPMENT  AND  PROCEDURE 

The  series  of  crystals  were  rectangular  quartz  bars  four  times 
greater  in  length  than  width,  with  the  thickness  adjusted  for  the 
desired  resonant  frequency.  Suitable  electronic  equipment  was  built  to 
drive  the  unplated  crystals  and  record  the  decaying  mechanical  vibra- 
tions. Figure  1 is  a picture  of  the  apparatus,  emd  a block  diagram 
of  the  electronic  system  is  given  in  Figure  2.  The  equipment  consists 
of  three  main  parts. 

3. 1 Experimental  Equipment 

3.1.1  Pulse  Forming  Network.  A pulse  forming  network  was  built 
to  give  pulses  of  variable  length  and  separation.  The  network  inter- 
rupts the  continuous  oscillator  signal  with  the  aid  of  a relay,  and 
provides  20  milliseconds  of  pulse  with  an  interruption  of  100  milli- 
seconds, which  is  long  enough  to  record  the  complete  waveform  on  a 
storage  oscilloscjpe.  The  generator  is  connected  parallel  tc  the 
crystal  at  all  times.  The  signal  voltage  applied  to  the  crystal  is 
60  dB  above  its  closed  gate  level. 

The  network  consists  of  three  IC  chips  and  a Babcock  relay.  A 
closed  loop  CD  4001  AE  multivibrator  and  its  external  circuit  controls 
the  repetition  time  and  provides  a 5-volt  square  wave  triggering  source 
for  the  9601  monostable  vibrator,  whose  external  circuit  controls  the 
pulse  width  and  drives  the  relay.  A peripheral  driver  is  inserted  to 
increase  the  current  to  the  relay.  A five  pf  capacitor  is  inserted 
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! 
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Figure  1.  Pictori"<l  of  the  equipment. 


( 
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between  the  relay  and  the  crystal  to  reduce  decoupling  the  crystal 
whenever  power  is  absorbed  by  the  electrical  equipment. 

3.1.2  Ten^erature  Bath.  A twelve-gallon  pyrex  glass  container 
with  suitable  heaters,  stirrers,  a thermocouple  gauge,  and  cooling  coils 
was  built  to  control  the  temperature  of  the  Ostwald-Fenske  viscometers 
whenever  measurements  at  larger  liquid  volumes  were  needed.  The  central 
portion  of  the  bath  was  vacuum  sealed  amd  the  essential  parts  of  the 
shear  viscometer  were  inserted  and  used  to  determine  the  correction 
factor  (see  Section  4.1.3).  Measurements  were  also  made  of  the  effect 
of  the  unavoidable  presence  of  adsorbed  water  on  the  crystal  surface. 

The  loss  factor  readings  proved  to  be  repeatable  in  both  air  amd  vacuum. 
Therefore,  the  shear  viscometer  was  removed  from  the  bath  and  placed  on 
a separate  stamd  to  allow  easier  accessibility  during  measurement. 

3.1.3  Pulse  Amplifiers.  A tuned  pulse  amplifier  was  built 
as  an  amplitude  limiter  to  prevent  the  2-voit  driving  pulse  from 
saturating  the  20  mv  crystal  output,  since  the  same  arrangement  is  used 
for  driving  amd  receiving.  Figure  3 shoivs  the  circuit  of  the 
amplifier.  Two  silicon  diodes  in  the  input  stage  prevent  the  incident 
pulse  from  blanking  out  the  amplifier.  The  2N2^22  transistors  provide 
a gain  of  680  over  the  frequency  range  from  200  to  620  kHz.  There- 
after, the  gain  drops  at  a rate  of  12  dB  pei  cctave  largely  due  to  the 
integrated  circuits  in  the  pulse  forming  network  limiting  the  anf)lifier 
to  frequencies  below  2.0  MHz.  The  amplifiers  amplitude  response  was 

determined  by  plotting  — in  dB  versus  time,  where  A is  the 

A + 1 n 

n 

decay  amplitude  n milliseconds  after  the  start  of  the  decay.  The 
straight  line  response  was  indicative  of  a linear  amplitude  response 
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over  a wide  range  of  amplitudes.  In  addition,  the  decay  time  of  the 
amplifier  was  less  than  1/10  of  the  normally  measured  decay  times. 

A demodulating  detector  forms  the  last  stage  of  the  amplifier.  The 
low  impedance  of  the  detector  to  high  frequencies  permits  the  low 
frequency  envelope  to  be  displayed  on  the  oscilloscope  screen.  The  time 
constant  of  the  detector  was  chosen  so  that  the  row  frequency  envelope 
coincided  with  the  output  of  the  last  tr^msistor  stage. 

A differential  pulse  amplifier  was  built  and  used  to  extend  the 
measurement  frequency  to  10  MHz.  The  silicon  diodes  again  serve  as 
amplitude  limiters  to  prevent  overloading  due  to  the  driving  pulse. 

The  WA733  operational  amplifiers  give  a fixed  gain  of  800  for  the  two 
stage  amplifier.  Figure  4 shows  the  circuit. 

3.2  Experimental  Procedure 

3.2.1  Procedure  During  Measurement.  Crystal  specimens  were  made 

by  Valpey  Fisher  Company  to  specifications  required  in  this  study  (see 

Appendix  B) . The  crystal  was  prepared  for  measurement  by  three 

repeated  washings  with  ethyl  alcohol  and  acetone.  This  procedure 

insured  that  the  surface  was  free  of  any  contamination  that  could 

prevent  the  crystal  from  being  in  good  contact  with  the  liquid.  The 

crystals  were  excited  while  resting  on  a bottom  base  electrode,  with 

the  top  electrode-crystal  separation  approximately  5.0  mm.  A drop  of 

liquid  was  then  placed  on  the  surface  with  a hypodermic  syringe,  amd 

the  decay  time  recorded.  Since  the  deciy  was  exponential,  the 

amplitudes  could  be  compared  from  the  slope  at  any  two  points.  The 

difference  between  this  time  constant  (t  ) and  the  time  constant  when 

e 
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the  cr/stal  is  driven  in  air  (t°)  is  a measure  of  the  effect  of  the 

e 

2 

liquid  viscosity  and  yields  a chemge  of  loss  factor  An  = . 

e 

The  crystal  was  earlier  placed  into  a vacuum  vessel  to  eliminate 

sound  radiation  and  the  relative  energy  of  the  dilatational  component 

of  vibration  determined.  The  resonant  frequency  w was  measured  with  a 

GR  1163-A  coherent  decade  frequency  synthesizer  to  an  accuracy  of 

1/10  Hz.  The  liquid  surface  contact  area  a was  measured  with  a micro- 

L 

scope  with  an  accuracy  of  0.01  cm^.  The  crystal  mode  meiss  was  the  final 
value  needed  ir.  Equation  (2.27),  amd  was  determined  (see  Section  4.1.5) 
to  be  — ^ times  the  crystal  mass  M®. 

2/7 

As  an  example  of  a particular  measurement,  the  decay  time  in  air 
was  measured  as  (58  i 1.0)  ms  at  502  kHz,  giving  a loss  factor  of 
(1.09  - 0.02)  X 10  *.  Placing  S-20  oil,  of  contact  area  0.155  cm*  , on 
the  crystal  surface,  decreases  the  decay  time  to  (22  i 1.0)  ms  and 
increases  the  loss  factor  to  (2.87  i 0.14)  ^ 10"® _ The  effect  of  the 
oil  increased  the  loss  factor  by  An  = (1.78  ± 0.12;  ^ 10"®.  For  502 
kHz,  the  static  mass  of  4.35  gms  gave  a mode  mass  of  1.53  gms.  ■,  The 
correction  factor,  from  Tcibie  II,  was  1.63.  Thus,  the  viscosity  of  the 
oil  droplet  was  (39.4  - 5.4)  cp.  The  uncertainty  of  the  viscometer  was 
14%  in  this  particular  measurement.  In  most  cases,  the  accuracies  were 
better. 

3.2.2  Procedure  for  Red  Blood  Cell  Deformability  Studies.  Bovine 
blood  was  used  as  the  source  of  red  blood  ceils  in  this  study.  The 
blood  was  collected  immediately  after  the  aorta  was  cut,  and  one  part 
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of  anti-clotting  solution  consisting  of  13.4  gms  of  sodium  oxalate  in 
1.0  liter  of  water  was  added  to  nine  parts  of  whole  blood.  The  whole 
blood  was  then  centrifuged  at  4000  rpm  for  15  minutes  so  that  the 
heavier  red  cells  could  be  separated  from  the  lighter  plasma  with  a 
thin  layer  of  white  blood  cells  remaining  in  between.  Both  the  plasma 
and  leukocyte  layer  were  syphoned  off  and  discarded  while  the  red  cells 
were  placed  in  am  equal  volume  of  saline  solution  and  washed  twice.  The 
washed  red  cells  were  then  combined  with  a normal  saline  suspension,  amd 
the  hematocrit  adjusted  for  desired  cell  voltune.  This  was  especially 
important  since  the  absorption  of  shear  waves  cam  be  attributed  to 
proteins  in  the  j.ed  cells  and  hemoglobin  is  the  most  ad>undant  protein 
in  blood  (10) . 

Chemicals  were  used  to  change  the  deformability  of  the  erythro- 
cytes. Sodium  Oleate  changed  the  rheological  properties  of  the 
erythrocytes  through  changes  in  cell  shape,  resulting  in  changes  in  the 
surface  area  to  volume  ratio.  Acetaldehyde  xs  known  to  increase  the 
tension  of  the  outer  membrane  of  a red  cell  with  a consequent  reduction 
in  surface  area.  The  preparation  procedure  was  as  follows:  1)  Two 

grams  of  sodium  oleate  were  added  to  100  milliliters  of  a normal  saline 
solution,  and  5 ml  of  this  solution  was  added  to  500  ml  of  red  cells 
resuspended  in  saline.  2)  Ten  milliliters  of  acetaldehyde  were  added 
to  324  ml  of  normal  saline.  This  solution  was  then  added  to  666  ml  of 
similarly  packed  red  blood  cells  to  make  one  liter  of  1%  acetaldehyde 
suspension. 

Due  to  the  great  effect  of  hemolysis  on  deformability,  a time  limit 
of  8 hours  was  established,  during  which  all  measurements  were  performed 
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The  frequency  behavior  of  droplets  of  blood  containing  these 
resuspended  cells  was  measured  with  the  shear  viscometer  by  measuring 
the  decay  time.  The  results  were  compared  to  available  Brookfield 
data  (11)  at  the  physiologically  low  shear  rates  from  0ol5  to  73  1/sec. 
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CHAPIER  :V 
EXPERIMENTAL  RESULTS 

4.1  Basic  Research 

4.1.1  Frequency  Response  and  Quality  Fac^c.:  cl  the  Crystals.  The 
frequency  spectrum  was  examined  for  crystals  vibrating  in  air  at  their 
resonant  mode#  with  a typical  result  given  in  Figure  5 for  283  kHz. 
The  figure  shows  a clear,  single  response  of  the  amplitude  and  its  ac- 
companying sidebands.  Spurious  modes  are  absent  within  40  Hz  on  either 
side  of  the  fundamental.  This  is  a critical  result,  since  adjacent 
modes  have  their  own  velocity  distribution  and  dilataticnal  components 
The  troublesome  neighboring  modes  have  been  elimi’^ated  by  proper  dimen- 
sioning, as  discussed  in  Appendix  B. 

The  figure  also  shows  the  respcnis  of  the  cristal  with  a drop  of 
NBS  oil  S-20  ^static  viscosity,  38  ? cp'  cr  _ts  surrace.  The  oil 
decresses  the  resonant  freguenty  by  10  Hz  and  suppresses  the  sidebands 
completely . 

To  insure  that  contribtt..ons  oi  modet  whose  frequencii  is  different 
from  the  driving  force  may  oe  neglected,  s :iystal  vibrat.ng  in  air 
should  have  a *arge  quality  factor  of  resonance  Q.  The  quality  factor 
of  the  cry  stals  wa=  plotted  from  0,1  to  0.5  MHz  at  ihr*r  lur.damental 
frequencies  and  from  06  to  1 0 MHz  at  an  odd  harmonic  of  the  fundamen- 
tal frequency.  A typical  result  was  iO'  at  0,5  MHz 

4.1.2  Support  of  the  Crystai.  5u*tab.e  ciystai  s..pports  were 
needed  to  insure  -nrestra-ned  vibra.-^rr-  -i  a*l  shear  modes  and  prevent 
the  generation  of  a secondary  f-aqjer.cy  spectrum  Ar  un?*-afcle  crystal 
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mounting  system  could  also  cause  a lack  of  reproducibility  during  mea- 
surement. The  quality  of  a certain  suspension  can  be  ascertained  from 
a plot  of  the  amplitude  of  vibration  for  various  relative  changes  of 
the  frequency  Aoj/W.  The  results  for  razor  blades,  phonograph  needles, 
and  a plane  metal  base  are  shovm  in  Figure  6'  The  figure  shows  that 
the  best  support  arrangement  is  to  simply  rest  the  crystal  on  the  brass 
base  with  the  crystal  essentially  suspended  by  an  elastic  air  cushion o 
Furthermore,  the  microscopic  air  layer  eliminated  sound  radiation  from 
the  bottom  crystal  surface. 

4.1.3  Effect  of  the  Electrodes  on  the  Vibration  of  the  Crystal. 

The  lower  electrode  is  termed  by  the  metal  base  of  the  crystal  holder. 
The  upper  electrode  is  electrically  separated  from  the  bottom  base  by 
four  insulated  support  rods  and  can  be  adjusted  vertically  by  a micro- 
meter screw  accurate  to  .01  mm.  The  air  gap  between  the  crystal  working 
surface  and  the  top  electrode  acts  as  a low  loss,  rigidly  reflecting 
wave  guide. 

Standing  waves  can  be  expected  to  influence  the  quartz  vibration 
by  introducing  interferences  that  depend  on  the  distance  of  separation 
Figure  7 is  a plot  of  the  quality  ractoi  'eusus  electrode  separation, 
from  0.25  to  30  mm,  for  three  different  loading  conditions,  vacuum,  air 
and  water  droplets..  The  maxima  in  th*s  i.gure  are  caused  by  ...nterfer- 
ences  from  standing  waves.  The  difference  in  the  value  of  Q when  the 
crystal  was  very  close  i0.25  mm,  so  that  sound  radiation  was  eliminated) 
to  a position  very  far  away  U.e.,  30  nuri)  is  a rough  indication  of  the 


amount  of  sound  radiated- 
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A more  accurate  measurement  can  be  obtained  by  examining  the  change 
in  the  quality  factor  in  air  and  under  300  microns  of  vacuum  when  the 
top  electrode  is  very  far  away.  (a. g.,  30  mm).  The  correction  factors 
in  Equation  (2.27)  would  thus  be  .determined  foi:  all  crystals.  Table  II 
shows  that  the  results  for  0.5  MHz  was  A = 0.33  and  *1  + A‘)^  = 1.23. 
Crystals  numbered  2 and  4 were  rejected  because  the  dilatational  compo- 
nent relative  to  shear  was  too  great. 

4.1.4  Effect  of  Pulse  Duration  and  Temperature.  The  influence  of 
the  relay  on  the  decay  measurements  was  considered.  At  the  moment  the 
power  is  switched  off,  the  relay  couid  be  at  a point  of  maximum  veloc- 
ity or  maximur.  amplitude,  or  somewhere  in  between,  The  amplitude  decay 
-6t 

Ae  may  therefore  depend  on  the  phase  of  the  relay. 

The  switching  on  and  off  also  generate*  undes..red  transients  To 
see  whether  they  have  any  effect  on  the  decay  amplitude,  the  pulse  dura- 
tion was  varied.  Figure  8 illustrates  the  decay  curves  for  different 
pulse  duration  irrespective  of  the  phase  at  the  moment  the  power  is 
switched  off.  The  linear  response  cf  the  amplitude  indicates  that  tran- 
sients and  phase  cf  power  switch  off  do  not  influence  tne  results. 

The  -ntluence  of  increasing  the  temperature  on  the  viscosity 
measurements  was  examined  by  varying  the  time  of  measurement.  Figure 
9 shows  the  results  for  drops  of  water  at  0 5 MHz  The  decrease  in 
viscosity  at  200  and  500  seconds  after  the  drops  have  been  placed  on 
the  crystal  surface  is  due  to  the  increased  temperature  resulting  from 
the  large  number  of  pulses  striking  the  crystal  Both  drops  approach 
a thin  film  geometry  because  of  evaporetion  as  seen  from  rhe  maxima  in 
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Table  II.  Ratio  of  dilatational  component  to  shear  component,  and 
correction  factor  (i  + versus  frequency. 


Crystal  # 

Frequency  (kHz) 

V , 
sh 

(1  + 

1 

235 

0.52 

1.63 

2 

282 

1.25 

6.65 

3 

283 

0.57 

1.74 

4 

296 

1.44 

9.47 

5 

402 

0.33 

1.23 

6 

503 

0,32 

1.23 

1 

602 

0,55 

1.69 

1 

704 

0.51 

1.57 

1 

997 

0-32 

1.21 

the  viscosity  curves.  Care  was  taken  to  finish  all  measuiements  with- 
in 200  seconds. 

4.1.5  Effect  of  Velocity  Distribution  Over  the  Crystal  Surface 
It  is  a very  difficult  task  to  determine  the  transverse  velocity  of 
small  crystal  bars.-  since  the  displacement  amplitude  is  of  the  order 
of  24  A/1  volt  at  500  kHz,  Tr.sre  is  also  no  reason  to  believe  that  the 
crystal  vibrates  with  a un*form  velocity  Pre.'icus  authors  (12)  have 
treated  the  crystal  velocity  distrioution  as  a c -e- dimensional  problem, 
with  the  displacement  varying  hanr.or loally  along  the  y aXiS, 

^ C;  sin(ky-o;t);  giving  an  effective  crystal  ma=s  cf  M = 4 The 

V ^ 

particular  transverse  velocity  distr-but-or.  also  determines  the  loss 
factor  for  the  liquid  through  its  dependence  cn  the  inverse  square  of 
the  liquid  surface  vaiOC-ty  amplitude 
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ProbeU>ly  the  nK)st  accurate  method  of  determining  the  mode  mass  is 
to  place  small  drops  at  various  positions  on  the  quartz  surface  and 
determine  the  local  velocity  from  the  resulting  decay  time.  A position 
of  larger  shear  motion  is  indicated  by  a short  decay  time,  and  a minimum 
of  motion  is  detected  by  a longer  decay  time.  The  loaded  crystal  wouxd 
have  a diminished  an^plitude,  yet  retain  the  displacement  pattern  of  a 
free  surface. 

This  has  been  done  for  glycerol  drops  and  the  results  are  repre- 
sented In  Figure  10.  For  a vibrating  quartz  bar  the  velocity  is  a 
maximum  at  the  center  and  decreases  towards  the  boundaries  of  the  plate. 
This  means  that  some  parts  of  the  plate  surface  take  no  part  in  the 
motion  and  the  average  velocity  must  be  determined  from  the  square  root 
of  the  local  viscosity.  The  one-dimensional  mode  mass  is  therefore 
closer  to  than  ^ M®.  Similar  distributions  were  measured  for  all 

2v^  2 

four  crystals. 

The  variation  of  velocity  with  width  was  determined  by  mapping  the 
two-dimensional  surface  of  the  crystal-.  A sharply  pointed  cone  was 
held  fixed  while  the  lower  crystal  support  was  driven  in  the  horizontal 
plane  by  two  clock  motors  through  a suitable  gear  arrangement.  The 
cone  functioned  as  the  top  electrode  as  well  as  a probe,  simultaneously 
applying  piezoelectric  excitations  locally.  The  crystal  surface  was 
mapped  by  synchronizing  the  signal  cf  the  stationary  upper  electrode 
with  an  oscilloscope  sweep  and  the  decay  amplitude  measured  as  a func- 
tion of  position  on  the  two-dimensional  surface.  The  central  region 
was  again  a maximum  of  shear  motion,  with  a drop  off  near  the  outer 
width  boundaries  influenced  by  the  decreasing  number  cf  flux  lines 
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rsachlng  the  crystal  surface.  The  results  show  that  the  vibrational 
velocity  does  not  vary  much  with  the  width  since  the  width  was  so  rela- 
tively small  compared  to  the  length.  Extensive  measurements  by  Sykes  on 
Y-cut  quartz  plates  have  shown  that  the  velocity  distribution  is  prac- 
tically independent  of  the  width  (13) . 

4.1.6  Effect  of  Contact  Area  of  the  Liquid.  Previous  authors  (14) 
have  found  that  energy  losses  in  liquids  were  significantly  influenced 
by  the  liquid  contact  area  In  this  stud^,  extensive  measurements 
were  made  of  viscosity  as  a function  of  liquid  contact  area,  and  a pat- 
tern began  to  emerge.  Figures  11  to  14  show  the  results  for  drops 
of  water.  The  peaks  in  the  viscosity  of  water  were  due  to  compressional 
resonances  that  generate  a relative  motion  parallel  to  the  surface  of 
the  crystal,  and  result  in  additional  shear  losses.  The  peaks  were 
therefore  ignored,  and  the  viscosity  was  determined  from  the  minima  of 
the  curves  after  repeated  measurements  were  made  on  each  sample.  The 
results  for  glycerin  and  urethane  varnish  are  shown  in  Figures  15  and 
16. 


4.1.7  Effect  of  Surface  Tension.  In  was  possible  to  increase  the 
volume  of  the  liquxd  without  changing  the  contact  area  by  placing 
smaller  drops  on  top  of  each  other,  in  this  manner,  the  radius  of 
curvature  and  the  effect  of  surface  tension  was  increased.  Also,  the 
surface  tension  could  be  decreased  by  adding  drops  of  detergent  to 
100  ml  of  water.  No  significant  change  of  the  viscosity  could  be 
observed,  although  the  viscosity  versus  drop  area  data  became  scattered. 
It  is  possible  that  the  method  used  to  incrsase  the  volume  (by  placing 
drops  on  top  of  each  other)  was  not  sensiti\e  enough  to  isolate  surface 


Figure  14.  Shear  viscosity 
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tension  effects  in  the  presence  of  coB$>ressional  droplet  resonances.  In 
fact  (see  Section  4.2)  v4ien  the  viscosity  was  later  measured  at  different 
frequencies r it  seemed  that  the  increase  of  droplet  viscosity  around 
400  )cHz  was  due  to  surface  tension..  This  increase  in  the  viscosity  at 
higher  frequencies  was  evident  for  all  the  liquids  thau  were  measured. 
Klemm  (15)  has  shown  that  surface  tension  could  increase  the  viscosity 
losses  by  a factor  as  great  as  10,000.  It  would  therefore  not  be  sur- 
prising that  a two-fold  increase  in  the  viscosity  could  occur  and  con- 
sequently superimpose  a scattezing  on  all  the  liquids  measured.  A 
greater  clarification  could  only  be  brought  about  by  building  equipment 
specializing  in  surface  tension  effects. 

4.2  Measurements  of  the  Viscosity  of  Normal  Fluids 

4.2.1  Water.  It  seemed  reasonable  to  believe  that  a calibration 
liquid  was  needed  to  determine  the  surface  losses  of  each  crystal,  and 
then  refer  all  other  liquid  results  to  this  fluid.  Work  performed  ear- 
lier with  glycerin  was  rejected  when  the  viscosity  of  glycerin  was 
found  to  be  strongly  dependent  on  the  water  content  and  varied  gr«  tly 
over  a relatively  small  temperature  range,  w^ter  was  finally  chosen  as 
the  most  suitaole  calibration  fluid  and  the  results  from  0.235  to  1.0 
MHz  are  are  shown  in  Figure  17.  The  viscosity  of  droplets  of  water 
is  close  to  the  static  value  of  0.89  at  25 -'C  for  all  values  except  in 
the  range  from  0.4  to  0.5  MHz.  The  increase  in  the  droplet  viscosity 
is  greater  than  the  uncertainty  of  the  measurement  and  is  not  due  to 
the  molecular  properties  of  tracer  since  any  relaxaticn  phenomenon 
should  show  a decrease  with  frequency,  is  likely  that  the  increase 
in  droplet  vi^scos  .ty  is  due  to  surface  tension.  The  surface  tension 
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effect  was  evident  in  all  liquid  droplets  that  were  studied,  except 
possibly  the  bonding  agents.  It  was  therefore  not  possible  to  refer 
the  experimental  results  to  water  and  absolute  measurements  were  pre- 
ferred. 

4.2.2  NBS  Oils  S-20  and  S-600.  The  viscosity  of  NBS  oil  S-20 
was  measured  over  the  frequency  range  from  0.283  to  1.0  MHz.  The 
results  in  Figure  18  are  scattered  about  the  static  viscosity  of  38.7 
centipoise  presumably  due  to  surface  tension. 

NBS  oil  S-600  was  measured  over  the  same  frec[uency  range,  and  the 
results  are  shown  in  Figure  19.  For  frequencies  below  0.4  MHz,  the 
results  are  near  the  static  value  of  1425  cp  at  25-'C.  However,  starting 
at  0.4  MHz,  the  viscosity  decreases  dramatically  and  scatters  about  rhe 
lower  value  of  400  cp.  This  decrease  in  viscosity  may  be  due  to  molec- 
ular relaxation,  which  is  possible  for  these  high  viscosities. 

4.2.3  varnish  and  Cedar  Wood.  The  viscosity  of  Satin  XL-8 
urethane  varnish  was  also  measured  over  the  frequency  range  from  0.235 
to  1.0  MHz,  and  the  results  are  shown  in  Figure  20  Below  0.6  MHz, 
the  viscosity  scattered  about  the  constant  value  of  73  cp,  slightly 
greater  than  the  static  value  of  60  ± 5 cp.  At  1.0  MHz,  the  viscosity 
was  19.5  cp.  This  large  decrease  in  the  viscosity  may  be  due  to 
absorption  of  energy  in  the  long  molecular  cr.ains  of  the  polymer  base 
of  urethane  varnish,  and  indicate  the  presence  of  relaxation  phenomenon. 

The  results  for  cedar  wood  oil  are  given  in  Figure  21 . The 
visc.'ity  was  greater  than  the  static  value  by  nearly  17  cp  possibly  due 
to  surface  tension.  Table  III  lists  the  results  for  water,  S-20,  S-600, 


varnish  and  cedar  wood  oil 


Figure  18.  Shear  viscosity  of  NBS  oil  S-20,  from  0.283  to 
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Figure  19.  Shear  viscosity  of  NBS  oil  S-600,  from  0.283  to  1.0  MHz 


Figure  21.  Shear  viscosity  of  cedar  wood  oil,  from  0.283  to 


4.2.4  Carbon  Disulfide.-  Xylene,  etc.  Other  liquids  covering  a 
wide  range  of  viscosity  were  investigated  at  a single  frequency  and 
the  results  compeured  to  the  static  viscosity.  TeUsle  IV  lists  the 
results  for  liquids  from  carbon  disulfide  to  xylene  at  0.283  or  0.402 
MHz.  For  most  instances,  the  viscosity  compared  favorably  to  the 
static  value . 

4.3  Measurements  of  the  Viscosity  of  Bonding  Agents 

Table  V gives  the  viscosity  for  commonly  used  bonding  agents 
ranging  from  caster  oil  to  transducer-fluid  28.  The  table  represents 
data  at  0.283  and  0.402  MHz,  and  the  dynamic  viscosity  is  compared  to 
the  static  values.  Of  all  the  liquids  measured  in  this  study,  these 
bonding  agents  (with  the  exception  of  glycerin)  had  the  least  amount  of 
scatter  in  the  viscosity.  The  viscosity  was  less  than  the  static 
values  for  all  the  agents  measured. 

4.4  Measurements  of  the  Viscosity  of  Body  Fluids 

The  viscosity  of  the  body  fluids  of  man  and  the  domestic  animals 
are  largely  unknown  due  to  the  small  volumes  of  most  fluids  that  are 
available  for  measurement.  With  this  in  mind,  bile  fluid  and  synovial 
fluid  were  withdrawn  from  bovines  and  fluid  was  extracted  from  the  eye 
of  a sheep  (all  three  fluids  were  obtained  from  the  Meat  Laboratory  at 
The  Pennsylvania  State  University)  and  their  viscosity  measured  at 
frequencies  from  0.283  to  0.502  MHz. 

Bile  is  a greenish,  viscid  fluid  that  is  secreted  by  the  liver 
to  neutralize  the  hydrochloric  acid  from  the  stomach-  The  synovial 
fluid,  as  mentioned  in  Chapter  I,  lubricates  the  joints.  The  results 
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Table  IV.  Shear  viscosity  of  numerous  liquids. 


Liquid 

Static  viscosity 
(cp) 

Frequency 

(kHz) 

Shear  viscosity 
(cp) 

Faint-thinner 

3.2 

283 

2.5 

Mineral  oil 

10.4 

283 

4.3 

Carbon-disulfide 

0.355 

283 

0.6 

Nitro-benzene 

2.03 

283 

0.8 

Xylene 

0.62 

283 

1.4 

Detergent  (Joy) 

20.1 

402 

12.3 

Table  V.  Shear  viscosity  of  numerous  bonding  agents. 


Liquid 

Static  viscosity 
(cp) 

Frequency 

(kHz) 

Shear  viscosity 
(cp) 

Glycerol 

622.0 

283 

574.0 

Caster  oil 

670.0 

283 

239.3 

Dow  Corning  200 

35.0 

402 

8.8 

Transducer  fluid  28 

100.0 

402 

73.4 

Transducer  fluid  X 

100.0  ± 5 

402 

71.3 
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for  bile  scattered  greatly  and  may  have  been  caured  by  sedimentation  of 
its  main  constituents,  producing  changes  in  the  viscosity.  The  viscos- 
ity of  joint  fluid  was  quite  regular,  as  were  the  results  for  the  eye 
fluid  which  was  primarily  made  up  of  water. 

The  viscosity  of  uncoagulated  whole  blood  taken  from  a finger 
prick  was  measured  at  two  frequencies,  and  the  results  shown  in  Table 
VI.  The  results  for  the  finger  prick  (consisting  of  whole  blood 
containing  fresh  red  cells,  leukocytes  and  platelets,  all  suspended  in 
plasma)  were  three  times  greater  than  the  viscosities  of  red  cells 
suspended  in  saline  that  had  undergone  metabolic  changes  during  storage. 

The  data  for  all  the  body  f.'uids  measured  are  shown  in  Table  VI. 
The  results  show  that  the  shear  viscometer  is  well  applicable  to  the 
measurement  of  liquids  of  special  interest  to  researchers  in  the 
medical  fields. 

4.5  Measurements  of  the  Effect  of  Chemical  Deformation  of  Red  Blood 

Cells  on  the  Viscosity  of  Blcod 

Blood  cells  are  presently  being  studied  by  a number  of  investiga- 
tors and  by  many  different  methods  (16) . The  easy  availability  of 
bovine  blood  cells  offered  a unique  opportunity  for  measuring  the 
response  of  the  viscometer  to  changes  in  the  rheological  properties  of 
the  red  cells  Specifically,  using  bovine  blood  as  a source  of  red 
cells,  the  deformability  was  changed  by  hematological  methods  and  the 
viscosity  measured. 

4.5.1  Effect  of  Hematocrit  on  Viscosity-  First,  the  effect  of 
red  cell  volume  (or  hematocrit)  on  the  viscosity  was  measured  for  two 


different  suspensions.  The  suspensions  were  normal  saline  and  normal 
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Table  VI.  Shear  viscosity  of  four  body  fluids  versus  frequency. 


Liquid 

Shear  viscosity  (cp)  at  the  following  frequencies  (kHz) 

283 

402 

503 

Bile 

2.4 

4.7 

11.0 

Eye 

1.8 

1.3 

Finger  prick 

15.4 

16.5 

Joint  fluid 

2.0 

1.8 

1.5 

saline  with  20  mg/100  ml  of  sodium  oleate  added.  As  seen  in  Table  VII, 
the  80%  cell  packing  has  a greater  viscosity  than  the  40%,  and  this  was 
true  for  both  suspensions.  Figure  22  shows  this  more  clearly.  As 
seen  in  Table  VII,  the  shear  viscosity  was  compared^le  to  the  measure- 
ments made  with  the  Brookfield  viscometer. 

4.5.2  Effect  of  Sodium  Oleate  on  Viscosity.  Sodium  oleate  was 
used  to  change  the  shape  of  the  red  blood  cells  and  the  effect  on  the 
viscosity  of  a suspension  containing  these  cells  was  measured  with  the 
shear  and  Brookfield  viscometers  (see  Table  Vlll)  . The  viscosity,  mea- 
sured with  both  instruments,  increased  when  20  mg/100  ml  of  sodium 
oleate  was  added  to  red  cells  suspended  in  normal  saline.  This  may  be 
seen  in  Figure  23. 

4.5.3  Effect  of  Acetaldehyde  on  Viscosity.  Acetaldehyde  offers 
a method  of  changing  red  blood  cell  deformability  without  introducing 
geometric  and  internal  viscosity  changes.  Therefore,  acetaldehyde  was 
added  to  bovine  erythrocytes  to  stiffer.  their  outer  membrane,  and  the 
effect  on  viscosity  was  measured  with  both  the  shear  and  Brookfield 
viscometers.  The  addition  of  1 ml/100  ml  of  acetaldehyde  increased  the 
viscosity  as  measured  with  both  instrxanents  (see  Table  IX) , The 
results  are  also  shown  in  Figure  24. 
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Figure  24.  Shear  viscosity  of  bovine  erythrocytes  in  normal  saline,  with  added  lml/100  ml  of 
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TaUsle  VII.  Shear  and  Brookfield  viscometrlc  measurei.'.ents  with  added 


20  mg/100  ml  of 

sodium  oleate 

versus  hematocrit. 

Shear  viscometer 

Shear  viscosity  measurements  (cp)  at 
the  following  frequencies  in  kHz 

Suspension 

H(%) 

215 

283 

402 

503 

Normal  saline 

40 

3.1 

2.9 

2.7 

2.7 

Normal  saline 

80 

3.7 

3.8 

3.3 

4.3 

Normal  saline,  with 
20  mg/100  ml  oleate 

40 

5.8 

4.1 

4.2 

4.3 

Normal  saline,  with 
20  mg/100  ml  oleate 

80 

5.8 

5.8 

8.0 

Brookfield  viscometer 

Viscosity  measurements 
following  strain  rates 

(cp)  at  the 
in  sec"^ 

7.3 

14.7 

36.7 

73.4 

Normal  saline 

40 

4.50 

4.10 

3.74 

3.59 
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Taible  VIII.  Shear  and  Brookfield  viscometric  measurements  of  bovine 
red  cells  resuspended  in  normal  saline ^ with  added 
20  mg/100  ml  of  sodium  oleate,  at  an  elapsed  time  of 
four  hours. 


shear  viscometer 

Shear  viscosity  measurements 
the  following  frequencies 

(cp)  at 
in  kHz 

Suspension 

H(%) 

215 

283 

402 

503 

Normal  saline 

80 

3.7 

3.8 

3.3 

4.3 

Normal  saline,  plus 
20  mg/100  ml  sodium 
oleate 

80 

5.8 

5.8 

8.0 

Brookfield 

viscometer 

Viscosity  measurements  (cp)  at  the 
following  strain  rates  in  sec~^ 

7,3 

14.7 

36.7 

73.4 

Normal  saline 

40 

5.00 

4.33 

3.91 

3.73 

Normal  saline,  plus 
20  mg/100  ml  of  sodium 

oleate 

40 

8.00 

6.69 

5.54 

5.00 
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Tcd;le  IX.  Shear  and  Brookfield  viscometric  measurements  for  bovine 
erythrocytes  resuspended  in  normal  saline,  with  added 
1 ml/100  ml  of  acetaldehyde. 


Shear 

viscometer 

Shear  viscosity  measurements  (cp)  at 
the  following  frequencies  in  kHz 

Suspension 

H(%) 

215  283 

402 

503 

Normal  saline 

40 

1.8  2.9 

2.7 

2.7 

Normal  saline,  plus 
1%  acetaldehyde 

53 

5c0 

5.6 

5.7 

Brookfield  viscometer 

Viscosity  measurements 
following  strain  rates 

(op)  at  the 
in  sec’^ 

7.3  14.7 

36.7 

73.4 

Normal  saline 

40 

4.70  4.38 

3.96 

3.83 

Normal  saline,  plus 
1%  acetaldehyde 

47 

6.00  5.65 

5.66 

5.90 
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CHAPTER  V 

SUMMARY  AND  CONCLUSION 

5.1  Suromary 

The  problem  of  this  investigation  was  to  develop  an  inexpensive/ 

precision  viscometer  capable  of  measuring  small  liquid  volvuneS/  such  as 

drops.  The  approach  was  to  utilize  the  short  effective  path  length  of 

shear  waves  in  liquids.  The  viscometer  consists  of  an  unplated  crystal 

bar  excited  by  a 50  ms  pulse,  and  suitcOsle  electronics  to  observe  and 

record  the  decaying  vibrations  after  the  power  is  shut  off.  There  are 

two  measured  time  constants  t defined  as  the  time  it  takes  for  the 

e 

waveform  to  reach  1/e  of  its  initial  amplitude.  The  first  decay  time 
is  for  the  crystal  driven  in  air.  The  second  decay  time  is  measured 
when  the  crystal  decays  with  a liquid  on  its  surface.  The  difference 
between  these  decay  times  is  a measure  of  the  effect  of  the  liquid  and 
gives  the  system  loss  factor. ' 

We  consider  the  case,  through  proper  experimental  arrangement, 
where  the  greatest  loss  is  the  shear  energy  loss  in  the  liquid.  All 
other  losses  are  negligible,  including  the  crystal  mounting  loss. 

A crystal  orientation  analysis  was  made  to  insure  a dominant 
thickness-shear  mode.  After  the  proper  crystal  cut  was  selected,  a 
dimensional  analysis  was  necessary  to  reduce  coupling  between  the  shear 
and  flexure  modes - 

There  was  no  reason  to  believe  that  the  crystal  surface  vibrates 
with  a uniform  velocity.  The  surface  velocity  distribution  was  found 
by  placing  drops  of  glycerin  at  various  positions  on  the  working 
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surface.  The  velocity  was  a maximum  at  the  center  and  decreased 
towards  the  bovindaries  of  the  plate.  The  crystal  mode  mass  was  deter- 
mined from  these  measurements  to  be  — ^ M®  , where  M®  is  the  crystal 
static  mass. 

Extensive  measurements  were  made  of  the  viscosity  as  a function  of 
liquid  contact  area,  since  the  energy  losses  were  influenced  by  the 
liquid  configuration.  Peaks  in  the  shear  viscosity  were  due  to  com- 
pressional  resonances  that  generate  a relative  motion  parallel  to  the 
surface  of  the  crystal.  These  peaks  were  ignored.  The  lower  parts  of 
the  viscosity  curves  were  averaged  ?.nd  the  average  constituted  the 
measured  viscosity  at  that  frequency. 

The  shear  viscometer  was  applied  to  a variety  of  liquids  ranging 
from  benzene  to  xylene.  The  measurements  for  a droplet  of  water  were 
close  to  the  static  value  of  0.89  cp  for  all  values  from  100  to  1000  kHz, 
except  between  400  to  500  kHz  where  the  viscosity  increased.  The 
ir'-rease  in  droplet  viscosity  was  greater  than  the  uncertainty  of  the 
measurement  and  was  not  due  to  the  molecular  properties  of  water  where 
a relaxation  phenomenon  would  decrease  with  frequency.  It  is  likely 
that  the  increase  was  due  to  surface  tension,  which  produced  a scatter- 
ing in  the  data  of  all  the  liquids  measured  in  this  study.  Unfortu- 
nately, the  method  of  measuring  the  effect  of  surface  tension  (by 
placing  drops  on  top  of  each  other,  thereby  increasing  the  radius  of 
curvature  while  keeping  the  surface  contact  area  fixed)  was  not  sensi- 
tive enough  to  isolate  surface  tension  in  the  presence  of  compressional 
droplet  resonances. 
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A National  Bureau  of  Standards  oil  S-600,  was  measured  from  283  to 
1000  kHz.  The  low  frequency  results  were  near  the  static  value  of 
1425  cp.  However,  beginning  at  402  kHz,  the  viscosity  decreased 
dramatically.  This  decrease  could  have  been  due  to  molecular  relaxation 
which  is  possible  at  these  high  frequencies  for  high  viscosity  oils. 

The  Satin  XL-8  Urethane  Varnish  viscosity  value  at  996  kHz  of 
19.5  cp  was  also  well  below  the  static  value  of  60  ± 5 cp.  The  decrease 
could  be  due  to  the  edjsorption  of  energy  in  the  long  molecular  chains  of 
the  polymer  base. 

A few  commonly  used  acoustic  bonding  agents  were  excunined,  ranging 
from  Castor  Oil  to  Dow-Corning  200.  The  data  for  these  bonding  agents 
was  the  most  repeatable  of  all  the  liquids  examined  in  this  study. 

The  shear  viscometer  orovided  a new  tool  for  examining  liquids  of 
special  interest  to  researchers  in  the  medical  fields.  Knowledge  of 
the  body  fluids  of  man  is  meager,  largely  due  to  the  small  volumes 
available  for  measurement.  For  instance,  the  viscosity  of  joint  fluid 
was  measured  and  the  result  varied  from  1.5  to  2.0  cp.  The  results  for 
uncoagulated  whole  blood  from  a finger  prick  were  three  times  greater 
than  blood  that  had  undergone  metabolic  changes  during  storage. 

The  applicability  of  the  shear  viscometer  to  medical  research  was 
examined  by  investigating  the  response  of  the  viscometer  to  rheological 
changes  of  bovine  red  blood  cells.  Specifically,  the  deformability  of 
red  cells  was  regulated  by  standard  hematological  methods  using  sodium 
oleate,  acetaldehyde  and  hematocrit,  and  the  viscosity  measured  with 
the  high  frequency  shear  viscometer  and  the  low  strain  rate  Brookfield 
viscometer.  Sodium  Oleate  changed  the  shape  of  the  red  cells  and 
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acetaldehyde  stiffened  the  outer  mendsrane.  The  results  show  that  the 
> dynamic  viscosity  increases  with  the  addition  of  sodium  oleater 

acetaldehyde,  or  hematocrit.  The  Brookfield  data  was  comparable. 

5.2  Conclusion 

Internal  friction  is  a consequence  of  molecular  and  atomic  relaxa- 
tion phenomenon.  For  fluids,  only  a limited  number  of  discrete  relaxa- 
tions are  superimposed  on  each  other.  The  viscosity  is  constant  at  low 
frequencies  as  long  as  the  sound  period  is  larger  than  the  releucation 
time  constant;  at  higher  frequencies,  the  viscosity  decreases  with  the 
frequency.  We  would,  therefore,  always  expect  the  viscosity  to  be  a 
maximum  for  static  flow  and  other  low  frequency  measurements  and 
^ decrease  only  at  the  higher  frequencies.  Therefore,  a peak  in  the 

viscosity  curve  indicates  an  unusual  behaviour  of  the  fluid.  Such 
peaks  have  been  observed  in  this  study. 

Most  measurements  in  this  study  lead  to  a value  of  viscosity  not 
very  different  from  the  static  viscosity,  and  some  (see  the  data  for 
NBS  S-600  oil)  indicate  clearly  a simple  relaxation.  However,  there 
are  some  instances  where  the  high  frequency  drop  viscosity  turns  out  to 
be  greater  than  the  static  value  (see  water,  for  instance) . This  could 
be  a surface  tension  effect.  It  is  still  unclear  how  this  effect  is 
coupled  with  the  sound  motion  and  with  the  radius  of  curvature  of  the 
drop.  It  is  possible  that  the  method  used  to  increase  the  volume 
(keeping  surface  contact  area  constant)  is  not  sensitive  enough  to 
isolate  surface  tension  effects.  It  has  also  been  suggested  that, 

, since  the  measurements  were  made  at  room  conditions,  the  liquid  layer 

near  the  drop  surface  accumulates  impurities  and  can  consequently 
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■i 

s 


prevent  the  crystal  from  being  in  perfect  contact  with  the  li^ui/d. 
Special  tests  could  be  designed  to  test  whether  impurities  eifeht  the 


measurements . 


The  variations  in  the  viscosity  were  greater  than  werM  ejected, 
and  detailed  studies  will  be  required  to  clarify  the  effec/n  of  surface 
tension.  However,  the  goal  of  this  thesis  to  produce  a ustfulL  method 
of  investigating  drops  of  fluids  has  been  fully  attained.; f Tae  method 
is  very  simple  and  practical,  and  the  results  are  reprodutilfLe. 
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APPENDIX  A 

CRYSTALLOGRAPHIC  ANALYSIS  OF  THE  CRYSTAL  'IBRATOR 

The  following  considerations  prove  the  advantag  -s  of  the  AT  quartz 
crystal  cut  over  the  normally  used  Y-cut, 

An  infinite  quartz  crystal  has  many  modes  of  motion:  flexural/ 

face-shear,  dilatational,  and  thickness-shear.  The  latter  three  types 

of  motion  are  coupled  by  mutual  elastic  coefficients  s^^  whose  value 

has  to  be  reduced  to  eliminate  these  troublesome  modes.  Consider  a 

special  infinite  quartz  plate  known  as  a Y plate, cut  such  that  the 

larger  face  contains  the  XZ  plane  emd  has  the  Y axis  for  a surface 

normal.  The  total  strain  induced  is  the  sum  of  the  mechanical  strain 

s, .T,  and  the  electrical  strain  d ,E,  , where  the  stress-strain  rela- 
ij  1 m3  j 

tions  exhibit  trigonal  symmetry.  The  total  electric  displacement  is 

the  sum  of  the  polarization  e .S  and  c E . The  sum  of  the  two  rela- 

nx  n mn  n 

tions  are 

Sj  = s.jTj  + Sj-Tj  + SjjTj  - S^^T^  - d. .E.  , 

^3  ~ ®33'^3  ' 

S4  = SjijTj  - SjijT^  + s^^T^  + d^i^Ej  , 


S5  - S^^Tj  2Sj^Tj  - dj^E2  , 


- t,^E^ 


2 ' 


D3  " ^3^.. 


and 


(A  1) 


where  Sj,  Sj,  S3  are  the  longitudinal  strains  and  S^,  Sj,  Sj  the  shear 
strains.  The  compliance  coefficients  s^g  (i  = 1,2,... 5)  couple  other 
inodes  to  the  thickness  shear  motion.  Of  special  importance  are  the 
coefficients  s..  and  s.  that  couple  dilatational  modes  to  shear.  In 
Y-cut  crystals,  Sjg  * s^g  = 0 and  dilatational  motion  should  not  exist. 

Consider  our  case  where  an  electric  field  Ej  is  the  initial  source 
of  strain.  Only  if  exists  in  a particular  row  in  Equation  (A.l) 
will  the  strain  of  that  row  survive.  Equation  (A.l)  becomes 


“ ^44*^5  ■*'  ^®14'^6  " *^14®2 


®6  “ 2sj^Tj  + 2 (s, |~S) j)Tg  “ 2d| jE 


11  *»12''^6  ^“11*2 


(A. 2) 


Dj  * "®mSj  - ®11®6  ■*’  ^11^2  • 


where  the  strains  are  referred  to  the  strained  coordinates  by 

9lj^ 

S,  * ^ (face-shear) 


5 9X3  3*1 


9ilJ  9\j; 

~ 9xi  ^ 9x2 


(A.  3) 


(thickness-shear)  . 


The  electric  field  Ej  therefore  produces  two  strains,  Sj,  Sg.  The 
Sj  strain  (face-shear)  distorts  the  crystal  in  the  xz  plane.  The 
Sg  (thickness-shear)  distorts  the  crystal  in  the  xy  plane.  Rewriting 
Equation  (A. 2)  in  terms  of  the  strains,  the  electric  coefficients  c^_ 
and  the  piezoelectric  constants  e^j. 


Ts  = 

''6  5^5 

CsgSg  + e-4E2 

Te  = 

®i 

(A  4) 

II 

CM 

0 

-6:455 

®ii56 
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Equation  (A. 4)  is  the  in^ortant  coupling  equation  between  face-shear 
> and  thickness  shear  The  strain  Sg  produces  the  Sg  strain,  the  elastic 

constant  Cgg  being  a measure  of  the  coupling, 

lb 

One  method  of  reducing  the  coupling  is  to  decrease,  by  a change 
in  the  orientation,  the  magnitude  of  Cggo  The  value  of  Cgg  depends  on 
the  rotation  of  the  crystal  cut  about  the  x axis  (that  is,  the  angle  9 
between  the  z and  the  rotated  z axis) . The  cut  of  0 * 31®  has  been 
called  the  Y cut.  The  cut  of  9 =*  35®,  knovrn  as  the  AT  cut,  has  the 
low  value  of  Cjg  * 60  ^ 10“^**  c.g.s.,  a zero  temperature  coefficient 
and  a relatively  active  electric  modulus  dj^  connecting  Ej  to  the  Sg 
strain  (17) . The  elimination  from  temperature  influence  is  especially 
important  when  using  shear  waves.  Crystals  from  the  AT  family  were 
chosen  for  this  study,  since  these  crystals  have  fewer  subsidieury 
1^  resonances . 

The  AT-cut  plate  also  gives  a maximum  of  energy  output  frcxn  the 
crystal  face.  The  place  is  cut  such  that  the  transverse  vibration 
has  its  motion  in  a direction  of  sub-minimum  elastic  modulus,  giving 
a maximum  amplitude  of  vibration . 
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APPENDIX  B 

DIMENSICftlAL  ANALYSIS  OF  THE  AT-QUAKTZ  PLATE 


When  dealing  with  plates  c£  finite  dimensions,  flexure  and  face 
shear  arise  from  the  presence  of  the  crystal  boundaries.  The  higher 
order  harmonics  of  flexure  are  especially  troublesome,  since  their 
equation  of  motion  is  identical  to  thickness  shear. 

It  becomes  necessary  to  select  dimensions  of  the  crystal  in  such 
a manner  that  the  three  types  of  modes  have  different  frequencies  and 
their  harmonics  do  not  intersect.  The  motion  of  a crystal  is  entirely 
determined  by  the  dimensions  of  the  plate.  The  coupling  of  the  un- 
wanted modes  can  therefore  be  avoided  by  choosing  the  proper  length  and 
width  such  that  tlie  unwanted  modes  do  not  intersect  the  experimental 
thickness-shear  frequency  curve  in  the  region  of  interest.  For  instance, 
the  regular  pattern  of  coupling  between  flexure  and  thickness-shear 
is  expressed  by  the  experimental  formula  (18) . 


,,  1338.4 

fv-  =■  — r~  , 

*f  X (mm)  *f 

where  x is  the  length  of  the  crystal  and  nxj  is  the  orde’^ 
for  n * 1,2, . . . 


(B.l) 

of  flexure 


The  coupling  to  harmonics  of  low  frequency  face  shear  depends 
on  both  the  length  and  the  widtn  dimensions,  Again,  experimental 
evidence  is  used.  The  frequency  of  coincidence  of  face  shear  and 
thickness-shear  in  AT  plates  in  expressed  by  two  separate  formulas 


(19)  , 
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and 


254.2 
X (cm)  *s 


kHz 


254.0 
z' (cm)”®s 


kHz 


(B.2) 


where  x * length  of  the  crystal, 

z'  ■ width  of  the  rotated  crj’stal, 
nj(^*  order  of  face  shear  along  x,  and 
nj5g*  order  of  face  shear  along  z' . 

Crystal  dimensions  were  predetermined  so  that  integer  flexure  and  face 
shear  were  not  present.  Table  X is  a list  of  the  crystals  bought 
from  Valpey  Fisher,  covering  the  ramge  from  100  kHz  to  5 MHz.  None  of 
the  crystals  have  an  integer  flexure  order.  Crystal  #7  does  have  a 
high  order  face  shear  (#13)  motion.  However,  this  order  is  too  high 
to  be  ‘oublesome. 
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Table  X.  Flexure  and  face-shear  modes  of  quartz  crystals 


Crystal (#) 

F(kHz) 

M(gms) 

Length 

x(mm) 

Width 

z(nin) 

Thickness 

y{mm) 

^Xf 

*s 

*s 

1 

106.5 

21.920 

45 

11 

17 

3.5 

1.9 

0.46 

2 

215.0 

10.060 

45 

11 

8 

7.2 

3.8 

0.93 

3 

283.8 

7.852 

45 

11 

6 

9.5 

5.3 

1.23 

4 

402.1 

5.45 

45 

11 

4 

13.5 

7.2 

1.74 

5 

502.9 

4.35 

45 

11 

3 

16.9 

8.9 

2.18 

6 

2000.0 

1.16 

45 

11 

0.9 

67.2 

35.2 

8.7 

7 

3000.0 

0.778 

45 

11 

0.6 

100.9 

53.2 

13.0 

8 

4000.8 

0.584 

45 

11 

0.45 

134.5 

70.8 

17.3 

9 

5000.5 

0.467 

45 

11 

0.4 

168.1 

88.5 

2'  .7 

) 
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APPENDIX  C 


DETERMINATION  OF  THE  CORRECTION  FACTOR 


Define  the  systexn  loss  factor  due  to  sound  radiation  to  be 


An  » -i-  Loss  due  to  Radiation  ■ 

u Total  Energy  of  the  System 

1 2 '^QPo®oVdil 


(C.l) 


2 


vrtiere  » surface  area  of  the  vibrator, 

PjCg  * characteristic  impedance  of  vibrator  in  air,  and 
Vdii  ■ maximum  dilatatJonal  velocity. 

1 

Approximating  the  mean  velocity  amplitude  by  ^ loss  factor 

becomes 


in-i?2£2^ 


(C.2) 


Since 


M - PgC^h 

* p c ~ , 
Q2  ' 


(C.3) 


where 


PqCq  » 42  gm/sec-CTi* 


pgCg  » 1.51  ^ 10  gm/sec-cm  , 


(C.4) 


therefore. 


1 + T2= 


1 _ 17.7  X 10 

An 


-6 


« 
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